A disordered RKKY lattice mean field theory for ferromagnetism in diluted magnetic 

semiconductors 
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We develop a lattice mean field theory for ferromagnetic ordering in diluted magnetic semiconduc- 
tors by taking into account the spatial fluctuations associated with random disorder in the magnetic 
impurity locations and the finite mean free path associated with low carrier mobilities. Assuming a 
carrier-mediated indirect RKKY exchange interaction among the magnetic impurities, we find sub- 
stantial deviation from the extensively used continuum Zener model Weiss mean-field predictions. 
Our theory allows accurate analytic predictions for T c , and provides simple explanations for a num- 
ber of observed anomalies including the non-Brillouin function magnetization curves, the suppressed 
low-temperature magnetization saturation, and the dependence of T c on conductivity. 
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Much of our current understanding of ferromagnetism 
in diluted magnetic semiconductors (DMS), most no- 
tably in the extensively studied molecular beam epi- 
taxy grown Mn-doped Gai-^Mn^As (with the Mn dop- 
ing level x « 0.01 — 0.1 system [1]), has been based on a 
very simple continuum Weiss mean field theory approx- 
imation [2] of the Zener model for the local (p-d) ex- 
change coupling between the impurity magnetic moment 
S = 5/2 e?-levels of Mn, and the itinerant carrier spin 
polarization, s — 3/2 holes of p-character in the valence 
band of GaAs. Spatial fluctuations associated with the 
random locations of Mn local moments are neglected in 
this continuum Weiss mean field theory (MFT). There is 
no theoretical rationale, except simplicity, underlying the 
neglect of the strong quenched spatial disorder due to the 
random magnetic impurity locations. Indeed, there have 
been several Monte Carlo simulations [3-7] attempting to 
include spatial disorder effects in the theory. These sim- 
ulations as well as a recently developed percolation the- 
ory [8,9] have explicitly demonstrated the manifest im- 
portance of quenched disorder in DMS ferromagnetism, 
at least for the localized insulating DMS systems. There 
are also strong experimental signatures for the impor- 
tant interplay between disorder and magnetism in DMS 
materials. 

In this Letter, we develop the first systematic theory 
for DMS ferromagnetism explicitly including spatial dis- 
order effects in the successful Zener-RKKY mean field 
model [2,9]. Our theory is explicitly constructed for 
the metallic DMS systems with itinerant carriers (since 
we assume the carrier-mediated effective Mn-Mn indirect 
magnetic exchange interaction to be of free carrier me- 
diated RKKY form), but the inclusion of a finite carrier 
mean free path in the theory (as we do in addition to the 
inclusion of Mn spatial disorder) allows us to make spe- 
cific predictions about the dependence of the magnetic 
behavior of the system on the carrier transport proper- 
ties. We also include a direct nearest-neighbor Mn-Mn 
antifcrromagnetic exchange interaction in our disordered 
RKKY-mean field theory, which takes on significance for 



larger Mn concentrations. An associated essential fea- 
ture of our model, which may be of crucial importance, 
is the fully discrete nature of our disordered mean field 
theory on a lattice. It is well-known that the discrete- 
ness associated with the specific lattice structure intro- 
duces unique features to the RKKY exchange interaction 
which are not caught in the corresponding continuum 
approximation. Our model, although conceptually sim- 
ple, is actually quite rich as it depends on five indepen- 
dent physical parameters of the DMS system: the carrier- 
induced Mn-Mn RKKY coupling (J ), the direct Mn-Mn 
(nearest-neighbor) antiferromagnetic exchange coupling 
(J ), the local moment density (rij), the free carrier 
density (n c ), and the carrier mean free path (ro). In 
principle, the randomness in the Mn locations on the lat- 
tice could also be parameterized, particularly if clustering 
(or other spatial correlations) of Mn impurities is impor- 
tant during GaMnAs growth. We neglect at this stage 
any such correlation in Mn spatial locations (since no 
independent experimental information on the nature of 
quenched disorder is available), assuming the Mn atoms 
to be uniformly randomly distributed in the zinc-blende 
GaAs lattice at the Ga substitutional sites. (It would be 
easy to incorporate a more detailed description of disor- 
der in the theory and to include Mn defects such as in- 
terstitials.) We also assume that the five parameters (i.e. 
Jo, J AF , rii, n c ,ro) of our theory are completely indepen- 
dent of each other, which is a theoretically and physically 
consistent approximation within our model since these 
parameters must be determined experimentally (or oth- 
erwise) and are independent parameters for our model. 

Our effective Hamiltonian describes the Mn-Mn mag- 
netic interaction between classical Heisenberg spins Si on 
a lattice: 
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where Si is the ith Mn local moment of spin 5/2; J, 
is the direct antifcrromagnetic exchange interaction be- 
tween nearest-neighbor Mn spins, i.e. J^ F = 0, unless 
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i and j are nearest neighbors. The first term in the ef- 
fective Hamiltonian, the carrier-mediated RKKY indi- 
rect exchange interaction responsible for producing DMS 
ferromagnetism, describes the magnetic interaction be- 
tween the Mn local moments induced by the extended 
(free) carrier spin polarization (these free carriers hap- 
pen to be the valence band or the impurity band holes 
in GaMnAs). This indirect Mn-Mn exchange interaction 
arises from [1,2,9] the local Zener coupling (or the so- 
called p-d hybridization) between the holes and the Mn 
d-levels, which then leads to the effective Mn-Mn RKKY 
interaction: 

jfj KKY = J r- 4 [sin(2fc F r) - 2k F r cos(2fc F r)] (2) 

where r = [Rj — Rj | is the spatial separation between the 

magnetically coupled Mn atoms; fc F oc nj^ is the Fermi 
wavevector corresponding to the carrier density n c , and 
the RKKY coupling strength Jo(> 0), taken as a param- 
eter in our theory, is related to the local Zener coupling 
J p d between the Mn local moments and the hole spins, 
J cx mJp d , where m is the hole effective mass. Note 
that for low carrier densities the RKKY exchange inter- 
action is mostly ferromagnetic except for large Mn-Mn 
separation (r > hp 1 ), and therefore, as long as n c <C rij, 
where rn is the active density of Mn local moments (i.e. 

— 1/3 

typical r oc n i ) , the frustration effects associated with 
the oscillatory nature of RKKY exchange interaction are 
unimportant in the problem; distinguishing the DMS sys- 
tems from random dilute metallic magnetic alloy spin 
glass systems (e.g. Cu-Mn), which are typically in the 
opposite limit of n c » rij. It is important to emphasize 
that a continuum Weiss mean-field theory (which neces- 
sarily neglects all spatial fluctuation effects by averaging 
over the Mn positions), blindly applied to the Mn-Mn 
RKKY interaction, as was already done [10] a long time 
ago (and revived recently [2] in the context of DMS sys- 
tems), always yields long-range ferromagnetic ordering 
of the Mn local moments for all values of n c and rii, 
with a mean-field ferromagnetic transition temperature 
rpMFT ^ j Qnin y 3 This is obviously incorrect for larger 
values of n c where ferromagnetism would eventually dis- 
appear [11] in a lattice model. 

We include the collisional broadening due to the finite 
carrier mean free path (rg) in the theory by incorporating 
an exponential suppression e~ r / r ° in Eq. (2), indicating 
that the RKKY interaction is cut off for r > ro due to re- 
sistive scattering effects. (This is important in view of the 
strong experimentally observed correlation between T c 
and conductivity typically found in DMS systems.) This 
form for the collisional broadening induced suppression of 
the RKKY interaction is theoretically well-justified [12] 
for studying ferromagnetic ordering of the local moments. 
We have also considered thermal effects in the RKKY 
interaction by generalizing Eq. (2) to the corresponding 
finite temperature formula, but typically T c <C T F for 



GaMnAs, where Tp (T c ) is the carrier Fermi temperature 
(Curie temperature), and therefore thermal broadening 
of the RKKY interaction is quantitatively and qualita- 
tively unimportant in our calculations in contrast to the 
collisional broadening effects included in the e~ r / r ° sup- 
pression. Note that an important novel feature of our the- 
ory is the "integrating-out" of the free carrier variables as 
we consider an effective (disordered) local moment mag- 
netic Hamiltonian defined by Eq. (1), with the free carrier 
information entering the theory only through the various 
physical parameters of the model, e.g. Jo,n c ,ro,m, etc. 
This simplification is well-justified in studying DMS fer- 
romagnetism since the long-range ferromagnetic ordering 
in the system arises entirely from the Mn local moments 
with the itinerant holes contributing little to ferromag- 
netism (by virtue of rij 3> n c and Smu > Shole)- 

The Hamiltonian of Eq. (1) can be rewritten as a gen- 
eralized Hciscnbcrg model for Mn (classical) spins on a 
disordered GaAs lattice H — J2ij Jij( r )^i ' Sj (where 
Jij = Jf- KKY + J{- F = J(r)) with the sum over i,j ex- 
tending over (random) magnetic impurity locations in 
the fee GaAs lattice. With no loss of generality we 
approximate the short-ranged antifcrromagnetic interac- 
tion jf- F by subsuming it in the definition of Jfj KKY or 
J(r) by multiplying it by a factor of 1/2 if the impuri- 
ties occupy nearest-neighbor positions in the fee lattice. 
This avoids introducing an additional unknown parame- 
ter J AF ', which could easily be introduced if such a need 
arises. 

We use a lattice mean field theory and consider each 
impurity spin Si to be immersed in an effective mag- 
netic field, fl<g = (M(S Z ))), where J t = £\ J tj is 
the sum of all couplings to the impurities surrounding 
the site i, and ((S z )) is the thermally and site-averaged 
polarization, and g is the g-factor corresponding to the 
impurity. Clearly, Jj depends on the impurity configura- 

(i) . . (i) 

tion and B c ^ is specific to the position Ri. Given B cf ^, 
the thermally averaged spin polarization (S?) is (fc F = 1) 

(Sf) = SBs (jjubB^/T^, where Bs is the usual mean 
field thermal Brillouin function. Although we calculate 
in mean field and use the average spin ((S z )), we retain 
the site dependence in Jj, thereby taking into account 
impurity disorder which subjects different Mn spins to 
different couplings depending on the local impurity con- 
figurations as determined by the random Mn locations on 
the GaAs lattice. (Sf ) is numerically calculated for a spe- 
cific impurity distribution. To determine ((S z )}, the po- 
larization averaged over all impurities, we integrate over 
all possible realizations of disorder, obtaining 

((S')) = JP(J)SB S (((S Z ))J/T)d.j, (3) 

where P(J) is the probability distribution of J. One then 
calculates ((S z )} self-consistently from Eq. (3). P(J) is 
determined a priori via Monte Carlo sampling. We as- 
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sume that the Mn impurities occupy Ga lattice sites with 
uniform probability as determined by its concentration x 
in Gai_ x Mn x As. 

In our lattice mean field theory the ferromagnetic tran- 
sition temperature T c is determined by the site-averaged 
exchange coupling and is given by 



i=l 



(4) 



where Ni and r; t are the numbers and distances of the 
ith nearest neighbors, respectively. We express and 
the remaining two length scales ro, (2k F )~ l in units of 
the fee unit cell length. The continuum limit is reached 
for r , (2k F y 1 < 1 with 

T cont^ (i4 07ra; /3) / r 2 J(r)dr (5) 
Jo 

= T c MFT [l-t a n- 1 (2k F r )/(2k F r )], (6) 

where T™ FT = (2807rJ a;/3)(37r 2 n c /2) 1/3 « J oXn l /3 is 
the continuum Weiss MFT value for T c which has been 
employed [2] extensively in the recent DMS literature 
(note that S = 5/2 here). The factor f{2k F r ) = [1 - 
tan -1 (2k F ro) / (2k F ro)] takes into account the exponen- 
tial cutoff imposed by the finite mean free path. We re- 
cover the continuum Weiss MFT result when 2k F ro 3> 1, 
i.e. in the strongly metallic regime. However, when 
ro becomes comparable to the length scale (2k F )~ 1 ; 

■ •, the RKKY interaction is effectively 
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suppressed and T c is substantially lower than T^t 1 even 
in the continuum approximation. To obtain an accurate 
formula for T c , it is necessary also to take into account 
the antiferromagnetic interaction between Mn impurities 
and to correct for the differences between the continuum 
approximation of Eq. (6) and the discrete lattice sum of 
Eq. (4). In the same way that an integral and a discrete 
approximation to that integral differ by a power series 
in the step size, the difference between the continuum 
and discrete formulae described above can be written as 
a series in k F . With these improvements, one finds as a 
reasonable large ro approximation for T c 

T r = T c MFT (l + a 2 n 2 J 3 + ai n 4 J 3 + a 6 n 2 c + • • •) 

- (1 - /?ir ->y 3 (l - /W)] + r„ - /WD 
f(2k F r ) - ^xJ(l/V2); (7) 
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x f(2k F r ) - — X J(1/V2); 
= -.733256, a 4 = 6.2594 x 1CT 2 , 
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a 6 = 2.89 x 1(T 3 , a s = 3.7 x 1CT 4 ; 
ft = .1633, & = .4284, f3 3 = .5584, 4 = .1176 (8) 

T c in Eq. (7) is correct to better than one part in 10 6 for 
large ro (as compared with our numerical calculations). 
Accuracy diminishes as r is made smaller, though the 
formula given in Eq. (7) is still correct to within 1% for 



ro = 1. One can also arrive at a formula to cover the 
small r regime. In this case, T c can be expressed as the 
sum of two terms 
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J(r)r 2 dr 



(9) 



where the first term in brackets is the exact nearest neigh- 
bor contribution while the second term contains the re- 
maining couplings calculated in the continuum limit, and 
the optimal choice for r a is 0.953. For r < 1 the rela- 
tive error for T c from Eq. (9) is less than 2%. Evidently, 
the first term (accounting for the effect of impurities on 
neighboring sites) dominates as r becomes comparable 
to the size of the GaAs unit cell. In this regime, the an- 
tiferromagnetic coupling removes half of T c (within our 
simplistic choice for J AF ) and thus plays a very impor- 
tant role for small mean free path. Obviously the result 
here depends strongly on the specific approximation one 
uses for J AF and would be stronger for stronger Mn-Mn 
AF coupling. 

In Fig. 1(a) we show our direct numerical calculation 
of T c as a function of the carrier density n c (for fixed 
Mn doping level x) for four different values of the carrier 
mean free path varying from strongly metallic (ro = 10) 
to "almost" insulating (r = .5). Clearly, for most choices 
of n c and ro, T^ FT is a poor approximation for T c in the 
disordered lattice system. We emphasize that the strong 
dependence of our calculated T c on the system conduc- 
tivity (through ro) is consistent with GaMnAs experi- 
mental results where increasing conductivity is found to 
enhance T c . In Fig 1(b) we show our calculated spon- 
taneous magnetization M(T) results, which depend on 
the full exchange distribution P{J). The M(T) profile 
(convex, concave, or linear) depends on whether the sys- 
tem is in the insulating (small ro, small n c ) or metallic 
(large ro, large n c ) regime. Concavity in M(T) is a signa- 
ture of an insulating system, while convex profiles appear 
deep in the metallic regime [9]. For intermediate impu- 
rity densities and mean free paths, it is possible to obtain 
a linear magnetization curve. Within the framework of 
our model, we are able to span both the localized and 
metallic regimes by having very small and large values 
of r , respectively, and we find both types of behavior in 
the magnetization profiles as can be seen in Fig. 1(b). 

There is a strong correlation between the degree of 
concavity in M(T) and the extent to which the coupling 
probability distribution P(J) has a strong weight near 
zero interaction strength and a multimodal profile. The 
multimodal P(J) in general leads to a concave M(T). A 
reasonable measure of the concavity is 7 = M (T)dT, 
or the difference in the slopes of M(T). The sign of 7 
indicates whether M(T) is convex (negative 7), concave 
(positive 7), or linear (if 7 w 0). The temperatures t\ 
and ti are chosen to capture an intermediate temperature 
range, neither very close to T c nor to zero. To go beyond 
the influence of the low temperature behavior, we set 
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FIG. 1. (a) 

rier mean free path ro = 10, 2, 0.9, and 0.5 Solid curves 
are plotted with the antiferromagnetic interaction taken into 
account, while dashed curves depict T c with antiferromag- 
netism suppressed, (b) Magnetization profiles correspond- 
ing to r = 10 (x = 0.05), r = 0.9 (x = 0.01), r = 0.5 
(x — 0.015); n c /n,i = 0.3 for all ro values 

t\ = .li c ; to avoid the critical region, we use t% — .7t c . 
This result is a rough criterion for concavity, since the 
sign of 7 indicates whether the profile is convex, concave, 
or linear. In Fig. 2 we give results for our "magnetization 
phase diagram" where the regimes of convex/concave 
magnetization behavior are depicted on the n c /rii—x two 
dimensional plots. Contour plots in Fig. 2 show in a 
concise way important trends; one sees very clearly the 
transition from some concavity in M(T) to convex M(T) 
behavior with increasing ro. Raising impurity density 
also tends to make the magnetization profile more con- 
vex. 

We have also calculated the saturation magnetization 
Mo = M(T — ► 0) using our theory, finding that consis- 
tent with experimental observations Mq could indeed be 
less than unity particularly for larger values of relative 
carrier density (n c /rii) and/or for more metallic systems 
(i.e. larger values of ro). The main suppression mecha- 
nisms are direct AF coupling between nearest-neighbor 
Mn-Mn interaction (increasing with J AF and x) and the 
oscillatory nature of the RKKY exchange coupling at 
large values of kpr (increasing with n c /n,-). 
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FIG. 2. Magnetization phase diagram: contour plots of the 
concavity parameter 7 as a function of Mn concentration x 
and ric/rii for two values of carrier mean free path (a) ro = .5 
and (b) ro = 2, where n c (rii) are the carrier (local moment) 
densities, and x is the Mn doping level. The 7 = contour 
divides the concave/convex phases in (a). 

To summarize, we have developed a disordered lattice 
mean field theory for DMS ferromagnetism which incor- 
porates spatial fluctuations associated with random lat- 
tice locations of the impurity moments, the finite carrier 
mean free path, and the Mn-Mn nearest-neighbor anti- 
ferromagnetic coupling. We calculate T c and magneti- 
zation curves for Gai-zMn^As ferromagnetic semicon- 
ductors, finding that all magnetic properties, including 
T c , depart significantly from the predictions of the con- 
tinuum Weiss mean-field theory. A particularly salient 
feature of our results is the strong theoretical correlation 
between T c and the metallicity of the system (i.e. ro) 
as observed experimentally. It should be feasible to in- 
corporate spin-orbit coupling and detailed valence band 
structure in the theory by generalizing our simple single 
valence band RKKY model. The most important essen- 
tial approximation of our model is the assumption of the 
RKKY form for the carrier-mediated indirect exchange 
interaction between the impurity local moments. Recent 
first principles band structure calculations [13] show that 
the effective DMS interaction between impurity moments 
is indeed of the RKKY form. There have also been recent 
numerical calculations explicitly establishing [14] the va- 
lidity of RKKY coupling in disordered semiconductors. 
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Our theory is computationally relatively efficient allow- 
ing us to obtain DMS magnetic properties as a function 
of five independent physical parameters (J , J AF , n c , n i} 
ro). We believe that the lattice mean field approximation 
is fairly sound by virtue of n, 3> n c in the DMS systems 
(and because of the long range nature of the ferromag- 
netic RKKY interaction). 
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